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ABSTRACT 



Context. The colour and metallicity gradients observed in spiral galaxies suggest that the mass-to-light ratio (M^,/L) of the stellar disc 
is a function of radius. This is indeed predicted by chemo-photometric models of galactic discs. 

Aims. We investigate the distribution of luminous and dark matter in spiral galaxies, taking into account the radial dependence of the 
stellar M^fL, which is usually assumed to be constant in studies of the mass structure. 

Methods. From earlier chemo-photometric models and in agreement with the observed radial profiles of galaxy colours, we derive 
the typical average M^/L profile of the stellar discs of spiral galaxies. We computed the corresponding variable mass-to-light (VML) 
stellar surface density profile and then the VML disc contribution to the circular velocity. We used the latter, combined with a well- 
studied dark matter velocity profile, to mass model co-added rotation curves. 

Results. By investigating rotation curves in the framework of VML stellar discs, we confirm the scenario obtained with the constant 
M*/L assumption: a dark matter halo with a shallow core, an inner baryon-dominated region, and a larger proportion of dark matter in 
smaller objects. However, the resulting size of the the dark halo core and of the inner baryon dominance region are somewhat smaller 
The stronger role that VML discs have in the innermost regions is important for constraining the galaxy mass structure in both A cold 
dark matter and MOND scenarios. 
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1. Introduction 

Within Newtonian gravity the mass model of spiral galaxies, in- 
cluding an (assumed spherical) dark halo, can be derived from 
high-quality rotation curves with little presence of non circular 
motions due to bars, spiral structure ot oval distortions; or from 
co-added rotation curves (e.g. Broeils 1992; Salucci & Burkert 
2000)13 A careful account of the distribution of the baryonic mat- 
ter is crucial. This includes the contribution of the stellar and HI 
discs, and of a significant central bulge in early type spiral galax- 
ies (not investigated here). 

The surface density of the stellar disc r*(^) is obtained from 
the surface brightness I(R), often fitted with an exponential pro- 
file of scalelength h (Freeman 1970) 



-R/h 



(la) 



by assuming that the stellar mass-to-light ratio (M*/L) is radially 
constant: 

lAR) = (MJL) m, (lb) 

where I{R) is either directly the observed surface brightness or 
its exponential fit. Equation [lb] corresponds to the simplest as- 
sumption according to which the stellar mass strictly follows 
the light. This is reasonable, especially when the surface bright- 
ness is measured in the red or infrared bands that best trace the 



' For the scope of this paper we use the term "rotation curve" also in 
the sense of "circular velocity as a function of radius". 



stellar mass distribution, in view of the approximate radial con- 
stancy of galaxy colours (e.g. Kent 1986; Verheijen 1997). On 
the other hand, various observational facts suggest that differ- 
ent regions in a spiral galaxy experience different evolutionary 
histories that lead to different colours and mass-to-light ratios. 
Then, it is worth dropping the assumption of a constant M*/L 
and investigating the effects of realistic radial variations in the 
M*/L ratio on the mass models. 

The main argument for this comes from the observed metal- 
licity gradients: in disc galaxies (including the Milky Way) the 
metallicity, traced by the HII regions and by the youngest stars, 
decreases outward with a typical gradient of -0.2 dex Z;^' in 
[O/H] (Garnett et al. 1997; van Zee et al. 1998; Shaver et al. 
1983; Smartt & Rolleston 1997; Gummersbach et al. 1998). The 
standard interpretation is that star formation and chemical en- 
richment proceed in the outer regions at a slower pace than in the 
inner regions, in an "inside-out" formation scenario (Matteucci 
& Frangois 1989; Sommer-Larsen 1991; Ryder & Dopita 1994; 
Portinari & Chiosi 1999 and references therein). This is sup- 
ported by the evidence that spiral discs are systematically bluer 
towards their outskirts and the scalelength of the light distribu- 
tion is typically shorter in red bands than in the blue (de Jong 
1996c; Bell & de Jong 2000; Jansen et al. 2000). Therefore, at 
a face value, stellar populations of younger age and of lower 
metallicity lie in the outer regions. All this points to a decline 
in M*/L with radius (Bell & de Jong 2001), whose extent and 
consequences are ripe for investigation. 
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The other main player in the mass structure of galaxies is 
the dark halo. This can be modelled by means of a theoret- 
ical approach that stems from the favoured A cold dark mat- 
ter (ACDM) cosmological scenario, very successful at account- 
ing for the large-scale structure of the Universe. For virialized 
objects, this predicts well-defined "cusped" profiles (Navarro, 
Frenk & White 1996) that, however, have been found to be at 
variance with observations (e.g. Moore 1994; Salucci, Walter & 
Borriello 2003; de Blok & Bosma 2002; Kuzio de Naray et al. 
2006; de Blok 2007; Spano et al. 2008). Alternatively, we resort 
to an empirical approach: one can easily find a form of (cored) 
DM density distribution that can account for all available kine- 
matic data (Gentile et al. 2004, 2005; Salucci et al. 2007). 

The present work aims to bridge the gap between the chemo- 
photometric studies, which support inside-out formation and ra- 
dial gradients in colours and M*/L, and the mass modelling of 
spiral galaxies where M*/L is traditionally assumed to be con- 
stant over the whole disc (with a few recent exceptions, see be- 
low). 

Let us briefly review the state of the art of the study of the 
kinematics of spiral galaxies with the aid of chemo-photometric 
information. The average M*/L ratio of a stellar disc, obtained 
from its global broadband colours by means of stellar population 
synthesis techniques, has proved to be a useful tool in dynamical 
studies of spiral galaxies (Ashman, Salucci & Persic 1993; Bell 
& de Jong 2001). More recently, Kranz, Slyz & Rix (2003) and 
Kassin, de Jong & Weiner (2006) for the first time have mass- 
modelled the rotation curve (RC) of a number of spiral galaxies, 
relaxing the assumption that M*/L(R)= constant and consider- 
ing the existence of (moderate) radial colour gradients. From 
these they derived the radial variation in M*/L by means of the 
(allegedly universal) relationship between the M*/L and colour 
of a mixed stellar population, 

< log >- a + b < (B - R) >, 

found by spectro-photometric models (Bell & de Jong 2001 ; Bell 
et al. 2003). This new treatment of the stellar disc surface den- 
sity, considered just as a refinement, has seemingly brought out 
no major new features in the mass models of spiral galaxies. 
Presently, colour-M*/L relations are being used in a variety of 
contexts where the mass profile of the stellar disc needs to be 
estimated (see e.g. Bakos, Trujillo & Pohlen 2008; Treuthardt, 
Salo & Buta 2008). 

However, this pioneering way of accounting for variable 
M*/L ratios is still subject to improvement. In fact, the above re- 
lation was derived from the global properties of model galaxies^ 
namely, the relation is supposed to hold among galaxies, relating 
their global colour and mass-to-light ratio. Yet, in the above- 
cited dynamical studies, the relation is instead applied within 
each object in order to estimate the local M*/L from the local 
colour {B - ^)(R). Moreover, it is sometimes applied beyond 
the range of colours for which it was established. Some caution 



^ The Bell & de Jong (2001) models are detailed multi-zone models 
accounting for the observed gradients in colour, metallicity, and age 
of galactic discs; however, the colour-M/L relations they present re- 
fer to their global galaxy models. Though similar relations have been 
shown to also hold radially, within individual discs (PST04) it is not a 
priori guaranteed that the same relations hold globally and locally, es- 
pecially for NIR bands where the age-metallicity degeneracy is lifted. 
This caveat applies in particular to the (much shallower) semi-empirical 
relations of Bell et al. (2003), which are the result of simple one-zone 
models calibrated on the global photometric properties of a large sam- 
ple of SDSS/2MASS galaxies. 



is therefore necessary when such colour-M*/L relations are ap- 
plied radially within individual galaxies: in the inside-out sce- 
nario, it is likely that radial variations in M*/L are far more "co- 
herent" than the variations related to global colour from galaxy 
to galaxy. 

In this paper we adopt a different approach to investigating 
the radial M*/L gradients and their consequences for the mass 
structure of spiral galaxies. We work out, for the first time in 
a self-consistent way, the theoretical radial gradients in M*/L of 
stellar discs by means of multi-zone chemo-photometric models. 
These follow the radial history of star formation and account for 
the observed metallicity and colour gradients. We obtain a typ- 
ical radial profile of variable mass-to-light ratio (VML) that, 
convolved with an exponential hght distribution, yields flie un- 
derlying stellar density profile. 

We derive the corresponding VML-disc contribution to the 
circular velocity and we mass-model the co-added rotation 
curves of Persic, Salucci & Stel (1996), that represent the av- 
erage RC of a spiral of a given luminosity well (see Salucci et al. 
2007 and references therein). Then, moving beyond the usual as- 
sumption of a constant mass-to-light ratio, we investigate crucial 
issues of the mass distribution of disc galaxies. These include the 
ubiquitous presence of DM in every object, the presence of a re- 
gion of "inner baryon dominance", the luminosity dependence 
of the dark-to-luminous mass fraction, and finally the inner den- 
sity profile of the dark matter halo. We also discuss VML discs 
in the MOND scenario. 

The paper is organised as follows. In Section 2 we derive the 
average typical M*/L profiles from chemo-photometric models 
of disc galaxies. In Section 3 we show how they are supported by 
the observed colour gradients. In Section 4 we determine the cor- 
responding VML-disc rotation curve; we then discuss the con- 
sequent disc/dark halo decomposition of rotation curves. Finally, 
in Section 5 we draw our conclusions. 



2. Chemo-photometric models and radial gradients 
in stellar mass-to-light ratio 

The standard "inside-out" scenario of disc galaxy evolution 
was originally developed to match a number of chemical prop- 
erties in the Solar Neighbourhood and in the Milky Way 
(e.g. Matteucci & Francois 1989; Chiappini et al. 1997, 2001; 
Portinari, Chiosi & Bressan 1998; Portinari & Chiosi 1999, 
2000; Boissier & Prantzos 1999). Subsequently, it was success- 
fully used to describe the chemical and photometric evolution 
of individual spiral galaxies (e.g. Molla, Ferrini & Diaz 1996, 
1997; Renda et al. 2005), as well as the whole spiral population 
(e.g. Boissier & Prantzos 2000, 2001; Prantzos & Boissier 2000; 
Bell & Bower 2000; Molla & Diaz 2005). In this section we dis- 
cuss how this scenario necessarily implies that stellar discs in 
spiral galaxies have radial gradients in colour and M*/L. 

We specifically consider the multi-zone models of Portinari 
et al. (2004; hereinafter PST04) that compute the detailed chem- 
ical and photometric profiles of the disc and the corresponding 
radial M*/L and colour profiles. Radial profiles were not explic- 
itly presented in PST04 (which was focused on the global M^,/L 
of disc galaxies), so we will do it in the present paper Let us 
anticipate ( ^2. 2b the main result: when tuned to observations, 
chemo-photometric models imply that the disc M*/L declines 
with radius R as 



-^(R) oc exp 
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with 0.5 < ,s < 1.1 and 0.14 </? < 0.8, and h the scalelength of 
the luminous profile. Any of the above cited works in the frame- 
work of the inside-out paradigm lead to similar M*/L profiles, 
though this fact is rarely mentioned explicitly (one exception is 
Fig. 7a of Boissier & Prantzos 1999), probably because the issue 
lies outside the direct scopes of chemical evolution studies. 

We briefly outline the main features of the code and the rel- 
ative models, but the reader uninterested in the details can di- 
rectly go to ^2. 2 1 where we present M*/L profiles derived from 
the PST04 models. In short; (1) the main observable underlying 
the derived mass-to-light profile is the metallicity gradient of spi- 
ral galaxies (-0.2 ±0.1 dex /i^' in [O/H], quite independent of 
mass and Hubble type); (2) the PST04 models correspondingly 
predict colour gradients that will be checked in Sect. 3. 

2.1. Main features of the models 

Here below we briefly summarise the main features of the 
chemo-photometric models of PST04, but further details can be 
found in the original paper. 

- The disc galaxy models possess a multi-zone radial struc- 
ture, divided into (typically 26) concentric annuli, which al- 
lows us to follow radial gradients in star formation (SF) his- 
tory. The local SF and chemical evolution history of each an- 
nulus are used to compute the radial profiles of metallicity, 
luminosity, and colours. 

- A galactic disc is assumed to form gradually via slow ac- 
cretion ("infall") of primordial gas, essentially cooling from 
the gaseous halo and settling into the disc. An exponentially 
decreasing infall rate is invoked to fit the metallicity distri- 
bution of long-lived stars in the Solar Neighbourhood (i.e. 
die "G-dwarf" problem; Lynden-Bell 1975; Tinsley 1980; 
Pagel 1997). This kind of infall is also observed in dynami- 
cal models of the formation of galactic discs (Larson 1976; 
Sommer-Larsen 1991; Burkert, Truran & Hensler 1992; 
Sommer-Larsen, Gotz & Portinari 2003). 

- The present-day surface density distribution of the baryonic 
matter (stars-Hgas) is forced to an exponential radial pro- 
file. With this choice, the present-day luminosity profile also 
approaches an exponential, as observations indicate, and a 
scalelength h can be confidently determined, at least out to 
3-4 scalelengths (see also Boissier & Prantzos 1999). 

- The code adopts a physically motivated star formation effi- 
ciency, which increases with the local surface mass density 
(Dopita & Ryder 1994; see Section 6 of PST04). This form 
of SF law is known to reproduce the radial chemical and pho- 
tometric profiles of galaxies (Portinari & Chiosi 1999; Bell 
& Bower 2000). 

- Once the stellar initial mass function (IMF) is fixed, the main 
parameters of the code are the SF efficiency v (in Gyr"') 
and the infall timescale t,„/. The models are then a two- 
parameter family. They were calibrated to reproduce the typ- 
ical average chemical properties (metallicity and metallicity 
gradients) of Milky Way-sized spiral galaxies, as follows. A 
grid of models with different infall timescales t,„/ is devel- 
oped, and for each Ti„f, the SF efficiency v is tuned to re- 
produce the typical metallicity of spiral galaxies with Mb ~ 
-21. Specifically, they reproduce 12H-log(0/H)=9.1 dex at 
a galactocentric radius R - 1 hs (see Sect. 6.2 of PST04). 
Once the absolute values of v and t,„/ are so calibrated, we 
consider the radial metallicity gradient that must match the 
observed -0.2 + 0.1 dex/i^' in [O/H] (Garnett et al. 1997; 
van Zee et al. 1998; Prantzos & Boissier 2000). This is of- 



ten obtained with no further tuning of the parameters. The 
adopted SF law in fact, with a dependence on surface density, 
can reproduce the observed metallicity profile. Only models 
with long infall timescales reach an agreement with observa- 
tions by introducing an additional radial increase in t,„/(7?) 
(TinfiR) oc R\ 0.5 <x< 2, see Table 2 of PST04). 
- We used six different IMFs: those of Salpeter (1955), Kroupa 
(1998), Kennicutt (1983), Larson (1998), and Chabrier 
(2001, 2002). The other IMFs are called "bottom-Hght" with 
respect to Salpeter, in that they store less mass in the low- 
mass end of the distribution (< 1 Mq) than the Salpeter IMF 
(extended down to 0.1 Mq). These IMFs are favoured by 
direct estimates of the stellar mass function in the Solar 
Neighbourhood, in the Galactic bulge, and in globular clus- 
ters and are indirectly supported by Tully-Fisher arguments 
(PST04). The upper mass limits on the IMFs were, in some 
models, tuned to reproduce the observed gas fractions of 
late-type disc galaxies, but such fine tuning has very little 
relevance for the aims of this paper. 

We thus obtained, for each of the IMFs considered, a series 
of models with (v, t,„/(/?)) calibrated to reproduce the observed 
metallicity of Milky Way-sized spiral galaxies, and the typical 
metallicity gradients (which are quite independent of mass and 
Hubble type). 

Increasing infall timescales t,„/'s "produce" objects with 
more extended SFHs (larger "birthrate parameters", the ratio be- 
tween present-day and past average SFR) and bluer colours, that 
represent galaxies of later Hubble types (Kennicutt, Tamblyn & 
Congdon 1994). 

Finally, we stress that, while the absolute value of metallicity 
is a function of galaxy mass, the metallicity gradient is quite 
independent of mass and Hubble type (Garnett et al. 1997; van 
Zee et al. 1998; Prantzos & Boissier 2000). Therefore, the radial 
trends obtained for the Milky Way-sized galaxies modelled by 
PST04 are actually valid in general. 

Barred galaxies may be an exception, as they often present 
much shallower metallicity gradients (e.g. Martin & Roy 1994), 
probably a result of bar-induced radial gas flows that wash out 
any pre-existing gradient. However, this phenomenon is prob- 
ably limited to strong bars; for instance, our Milky Way hosts 
a bar (or two; Lopez-Corredoira et al. 2007, and references 
therein) yet its metallicity gradient is in line with the one we 
model — and has actually inspired the whole inside-out scenario 
since the 80's. Besides, bar-induced gas flows wash out only the 
present-day metallicity gradient in the gas phase, while the un- 
derlying metallicity, colour, and M*/L gradients in the stellar 
component are much less affected, so our models should reason- 
ably describe barred galaxies, as long as their colour gradients 
are similar to those of other galaxies. 

Clearly, our present discussion does not apply in individ- 
ual objects where no significant metallicity and colour gradi- 
ents are observed (whether due to bar mixing or other reasons); 
however, this paper focuses on average trends and we show in 
Section |3] that the models yield reasonable predictions for the 
average colour gradients of disc galaxies. For the sake of com- 
pleteness, we mention that radial mixing of stellar populations 
may also affect the chemo-photometric evolution of disc galax- 
ies (Roskar et al. 2008; Schonrich & Binney 2009). While very 
promising for understanding the detailed connections between 
kinematics and chemistry in Galactic stars, models including 
these effects are still in their infancy. Whether they represent 
just a refinement over the underlying inside-out scenario or will 
drastically change the standard interpretation of the metallicity. 
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colour, and M*/L gradients of disc galaxies remains to be ex- 
plored, especially considering that such models involve many 
more free parameters than the standard ones. In this paper we 
necessarily limit our discussion to the dynamical consequences 
of the decade-long, well-established inside-out scenario. 

2.2. The radial profile of M^/L 

It is well known that the absolute value of M*/L depends on 
the SF history, stellar metallicity distribution function, and un- 
derlying IMF. Global information, such as absolute values of the 
metallicity and colour, is very valuable to frame the galaxy prop- 
erties and compare galaxies to each other — though is per se un- 
able to constrain the mass models with a precision comparable 
to kinematic analysis (see Salucci, Yegorova & Drory 2008). 

However, in this paper we are not interested in differences 
among different objects, but focus on the internal radial profiles 
in M*./L and their dynamical consequences compared to the sim- 
ple constant M*/L assumption. Therefore, we shall adopt a suit- 
able scaling that cancels out object-to-object offsets and high- 
lights the radial M*/L variation within the discs. 

We define the / band scalelength h; by fitting the model 
/ band luminosity profile with the exponential law of Eq. [la] 
in the range 0.5 - 3 hi. Introducing the notation M*/L^(7?) = 
we take the value of the M*/L at one scalelength, 
'^A,\ =M*/L^(/!/) as the zero-point normalization of the M*/L 
and we define the normalized M*/L radial profile as 



T 



(2) 



The normalization eliminates the zero-point offsets related to 
the specific IMF and star formation history of each model. The 
normalized M*/L profiles of the models show very similar be- 
haviour, with representative examples shown in Fig. [1] The pro- 
files range from "shallow" (close to linear out to ~ 3/i and shal- 
lower farther out; top panel) to "concave" (bottom panel). 

In any case, M*/L decreases with radius in all models, sim- 
ply as an effect of the age and metallicity gradients of the stellar 
populations that evolve in an inside-out scenario. Quantitatively, 
M*/L/ decreases by about a factor of 2 between 1 and 3^ scale- 
lengths, with a rise towards the centre for the concave cases. 
Though the radial gradients in M*/L are stronger in bluer bands, 
they are also non-negligible in the K band: M*/L/s: decreases by 
about 40% between 1 and 4/1/. As a result, the stellar surface 
density profile is generally more concentrated than the bright- 
ness profile, even in the NIR (see also Kranz et al. 2003). 

To proceed, we first consider that the / band offers a number 
of observational and theoretical advantages (see PST04), so that 
we choose it as the reference band — as in the kinematic anal- 
ysis of § 14.11 We find that the normalized M*/L/ profiles are all 
approximated well by 



Y 



= exp < -a{s) 



with 



a{s) = -(1.3 - i)^ -hO.13. 



(3a) 



(3b) 



In detail, shallow M*/L/ profiles have i ^ 1 (0.9-1.1), i.e. 
close to exponential, while concave profiles are well-fitted by 
s - 0.5 - 0.7; e.g., the example concave profile in the bottom 
panel of Fig. [1] corresponds to s = 0.5. The standard constant 
mass-to-light assumption is recovered for s = 0. 




shallow profile 



0.4 



2.5 



2 

R/h, 



ri.5 




^■^ _ concave profile 



12 3 4 

R/h, 

Fig. 1. M*/L profiles in different bands predicted by our chemo- 
photometric models, scaled to the respective values ?A R - hi 
(Eq.|2]i. An example of M*/L with shallow profile (model Kenn- 
H of PST04, top panel) and one with concave profile (model 
Chab-F, bottom panel) are shown. 



The analytical fit in Eq. 3 is accurate to 2-3% at /? > 0.5 h 
and to about 10% in the very central regions. For the typical 
"average" spiral galaxy the actual value of s will lie within the 
(reasonably narrow) range of values discussed above. In pass- 
ing, let us point out that it is up to future investigations to test 
whether Eq. 3 can be also representative of individual objects, 
and determine the corresponding values of s. 



3. Colour gradients 

It is well known that spectro-photometric models of spiral galax- 
ies predict tight relations between the M*/L and the colours 
of (simple and composite) stellar populations (Bell & de Jong 
2001; PST04). Henceforth, the observed colour gradients in spi- 
ral galaxies can be considered the empirical counterpart of the 
M*/L gradients. An obvious test for our model M*/L profiles 
is then to compare the corresponding colour profiles, predicted 
to get bluer outwards, to observations. To this purpose we con- 
sider the multi-band sample of 86 disc galaxies of de Jong & van 
der Kruit (1994). It is worth noting that the overall colour of a 
spiral galaxy tends to be redder for earlier morphological types 
(e.g. Roberts & Haynes 1994) because of their different aver- 
age disc SF histories, rather than to systematic variations in the 
bulge-to-disc ratios (Kennicutt et al. 1994; de Jong 1996bc). To 
circumvent this intrinsic variance in the global colour of differ- 
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R/h, 



Fig. 2. Colour difference profiles for the sample of disc galaxies 
of de Jong & van der Kruit (1994). T codes the galaxy mor- 
phological type: T=0-2 (Sa-Sab), T=3^ (Sb-Sbc), T=5 (Sc), 
T=6-10 (Scd and later types). Dots with error bars represent the 
average A(colour) and its 1 cr dispersion for the co-added pro- 
files. 

ent galaxies, and highlight the radial trends within each of them, 
we normalize the colour data by taking the colour measured at 1 
scalelength li as the zero-point and by shifting all the values as 

A(colour) - colour(^) - colour(/i). (4) 

This shift is consistent with the M*/L scaling of Eq. |2] that we 
apply to the models. 

We show the resulting "colour difference profiles" in Fig.|2l 
which is obtained from Fig. 2 of de Jong (1996c) after the zero- 
point shift by expressing the radial coordinate in disc scale- 
lengths (taken from de Jong 1996a). For the sake of example 
we show the run of B - K colour versus radius measured in R- 
band scalelengths, but similar trends hold in other bands (see 
also Fig.O. We ignore the data of the innermost half scalelength, 
which are heavily affected by the bulge component. Fig.|2]clearly 
highlights the existence of negative radial colour gradients, out 
to at least 3 scalelengths. For different morphological types the 
colour gradients are quite similar in the range 0.5 < R/h < 3, 
with Scd and later type galaxies showing a larger scatter. Beyond 
3h the colour profiles become shallower (and noisy). 

In spite of the considerable scatter (and of uncertainties in 
sky background subtraction at the outermost radii), there is no 
question that the colours on average become bluer outwards; 
B - K, for example, decreases by about 0.25 mag per disc scale- 
length. From this, we can claim that colour gradients emerge as 
a systematic, constitutive property of galactic discs. 

Then, in Fig. [3] we overlay the colour profiles of all morpho- 
logical types < r < 6 (neglecting later, more iiTegular types) 
and we compare their average values with the predictions of our 
chemo-photometric models. In detail, we plot various colour 
profiles in different bands and compare the data with the model 
predictions obtained by assuming a variety of IMFs. Each panel 



shows, for a given set of IMF models, two examples of model 
colour gradients, the shallower/steeper one corresponding to the 
shallow and concave M*/L profiles in ^2.21 

For all the adopted IMFs, a good agreement results between 
the colour profiles predicted by the models for an "average" spi- 
ral and the actual observations. Therefore, in the inside-out for- 
mation scenario of stellar discs, colour gradients emerge natu- 
rally and are remarkably similar to the observed ones. 

Notice that we neglect the innermost Q.5h region, where the 
comparison would be inappropriate, as data is affected by a bulge 
component, not included in our models. It is not our aim to study 
bulge-dominated discs or the bulge-disc interface, and since the 
innermost R < 0.5 h region encloses only 10-15% of the disc 
mass, the lack of a more exact treatment does not significantly 
bias our investigation on VML-discs. 

In conclusion, the observational properties of stellar discs 
indicate that their M*/L is likely a function of radius and that 
chemo-photometric models predict the average radial colour and 
M*/L profiles of spiral galaxies reasonably well. 

3.1. Dust effects 

In our models and in the comparison with the observed colour 
profiles, we have neglected the effects of dust. This is to a first 
approximation justified, since the observational sample consists 
of face-on galaxies for which dust has been found to play a mi- 
nor role in the colour gradients (de Jong 1996c). For an extended 
sample of low inclination galaxies. Bell & de Jong (2000) sim- 
ilarly conclude that colour profiles are driven by gradients in 
the age and metallicity of the stellar populations; dust induces 
a much smaller additional gradient. Besides, again to first ap- 
proximation, dust does not alter the colour-M*/L relation (the 
reddening-extinction dust vector running almost parallel to the 
relation, see Bell & de Jong 2001). Therefore the M*/L profiles 
coiTesponding to the observed colour profiles are also expected 
to be mainly driven by stellar population gradients rather than by 
dust effects. 

On the theoretical side, the chemo-photometric model of 
Boissier & Prantzos (1999) for the Milky Way also predicts that 
dust slightly enhances, but does not create colour gradients. As 
their models are quite similar to those of PST04, we take the re- 
sults from their Fig. 7 to illustrate the point. In Fig.|4]we show 
the Mi,fLv and the colour profiles of the Milky Way model of 
Boissier & Prantzos (1999): dust-free predictions and (face-on) 
extinction. The left panels reproduce their Fig. 7, with the ra- 
dial coordinate changed from R to R/hR, where hu ~ 3 kpc 
(estimated from the surface brightness and colour profiles in 
their Fig. 7). The 7?-band scalelength should be quite close to 
the /-band scalelength, which we refer to in our discussion of 
the M*/L gradients and rotation curves. The right panels show 
the same M*/Lv and colour profiles, but with the relative scal- 
ings adopted in Sects. l272l and[3] Clearly the absolute values of 
colours and M*/L are affected by dust (left panels), but in rel- 
ative terms, the effect on the colour gradients is small with re- 
spect to the amplitude of the gradient. The strongest effect is on 
optical-NIR colours like B - K, as noticed by de Jong (1996c). 
Comparison to Fig.[3]shows that the effect of dust on the colour 
gradients is much less than the scatter among individual galax- 
ies, and it does not significantly change the discussion in the 
previous section. As to M*/Lv, dust enhances the gradient but 
remains a second-order effect with respect to the main trend set 
by the inside-out scenario; and extinction affects the /-band (rel- 
evant for the remainder of our paper) much less than the V-band. 
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Fig. 3. Radial colour difference profiles in different bands, computed relative to the colour at R = h (Eq.|4]l. Thin dotted lines: data 
from de Jong & van der Kruit (1994) for < T < 6; dots with error bars: coiTesponding average A(colour) and 1 cr dispersion; 
thick solid and dashed lines: chemo-photometric models from PST04 (specifically, the models plotted are Krou-C,E; Salp-A,E; 
Chab-A,F; Kenn-H,J). 



In summary, dust certainly plays a non-negligible role when 
discussing the absolute values of colours, M*/L ratios, and age 
of the underlying stellar populations, but its role is secondary in 
terms of the gradients discussed in this paper. Dust effects cer- 
tainly deserve further attention, especially for inclined galaxies, 
but they will not change our main point. Radial M^/L gradients 
due to stellar populations must exist according to the inside-out 
scenario, and they significantly change the disc rotation curve 
with respect to the classic Freeman one (Sect. |4]i. From Fig. |4] 
we expect that dust at most enhances the trends we discuss. 



4. Variable mass-to-light rotation curve 

To establish the actual importance of the radial variations of the 
stellar M*/L, a complete chemo-photometric investigation of a 
large sample of individual objects would be needed, which is 
beyond our present scope. Here we aim to investigate the effects 
of the average stellar M*/L gradient of spiral galaxies (Eq. 3), 
in connection with their average rotation curve (RC), defined as 



the co-added curve of a suitably large sample of RC's of objects 
of fixed luminosity. 

From now on, we specifically consider the / band, which 
is the reference band for the set of rotation curves studied in 
§ 14. II and drop the subscript / for simplicity. From Eq. 3 for the 
M*/L profile we can build the average luminous component of 
the gravitational potential of spiral galaxies. We consider values 
of ,5=0.5, 0.7, and 1.0 (Fig.|5j, which cover the range of profiles 
obtained from the models. 

We assume that the Freeman profile (Eq. [Ta]i describes the 
light distribution. The (small) deviations we observe in the sur- 
face brightness of some spirals are important for individual mass 
modelling but irrelevant for the present purposes. As a conse- 
quence of the radial gradients in M*/L, the stellar surface density 
2*(/?) is somewhat more concentrated than the surface luminos- 
ity profile. For an exponential disc (Eg. [Tab, from Eq. 3 we get 



2*(7?) = 2i tx^{-a{s) 



(fr-'i-[(i) 



(5) 
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Fig. 4. Radial M*/Lv and colour profiles from the Milky Way model of Boissier & Prantzos (1999); thin lines: dust-free predictions; 
thick lines: including extinction. The right panels show the relative M*/Lv and colour profiles, scaled as discussed in § l2.2l and |3] 



where S] - l.-kih) - Ti lo/e is the surface density at 1 scale- 
length. For a given 2i, the surface density has a central value 
£*(0) - £ie"^', higher than that of the M,»./L=constant case by 
a factor e" 1.2-2 for s in the range 0.5-1.1), and decreases 
more rapidly outward than in the standard mass-follows-light 
assumption. In the case of 5 = 1, the surface density still has an 
exponential profile, but with scalelength 1 +a), i.e. about 20% 
shorter than that of the surface brightness profile. In this case 
we note that the main VML feature, a mass scalelength slightly 
less than the light scalelength, may have been so far masked by 
the intrinsic difficulties in deriving the surface brightness scale- 
length h well within the 20% precision level. 

In Fig.|5]we show three analytical M*/L profiles with i=0.5, 
0.7, and 1.0, and in Fig.|6]we show the corresponding cumula- 
tive disc masses compared to the constant M*/L (CML) case. 
The radial mass build-up is obviously more rapid and saturates 
more quickly in the VML case, with up to 70% of the total mass 
contained within about 2 h and up to 90% of the mass within 
Ropt =3.2 h. 

By solving the appropriate Poisson equation we compute nu- 
merically the self-equilibrium circular velocity of thin discs with 
a VML surface density as given in Eq. |5] In Fig. |7] we show 



the resulting rotation curves for the three representative VML 
discs of Figs. |5] and |6] compared to the standard CML rotation 
curve. The VML RC rises more steeply and peaks at smaller 
radii (1.5-1.9 h) than the Freeman CML rotation curve (2.15 h), 
and it falls off also more rapidly after the peak. The VML RCs 
approach the corresponding Keplerian fall-off (defined as the 
region where the logarithmic slope of the RC is close to -0.5) 
already at R = 3 - 3.5 /i, rather than at R > 4h as it occurs for 
the CML disc; fovR^Sh the VML RCs finally overlap with the 
CML RC. 

It is useful to provide, for the general VML RCs, an analyti- 
cal fit out to 3.5 h: 



■0.9+0.4.5 



(^2 + (i.35 + o.5i)2) 



2\l+0.3.! 



(6) 



with^ = R/h and n(s) = (1 -h(1. 35 +0.5i)^)'+°-^^" a normalization 
factor. The fit is precise to better than 2% between 0.3 < f < 3.5, 
and better than 10% down to ^ > 0.15. At larger radii, the VML 
RCs are described well by a Keplerian fall-off. The relation be- 
tween the total disc mass and the value of the circular velocity 
at R,,, - 2.15 h, the radius where the CML disc RC has its maxi- 
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Fig. 5. Three representative examples of analytical M*/L/ pro- Fig. 6. Cumulative disc masses corresponding to the M*/L/ pro- 



files derived from the chemo-photometric models (Eq. 3). 



files in Fig. |5] compared to the constant mass-to-light (CML) 
assumption. 



mum, reads as 

with 77 = (0.98, 1, 1.03) for the cases of s = (0.5,0.7, 1), and 
T] = 1.2 for the Freeman CML disc. 

4.1. The dark-luminous matter interplay 

The stellar surface density 2*(/?) in spiral galaxies is then only 
exponential to a first approximation; more in detail, it is an ex- 
ponential modulated by the M*/L gradient, as given by Eq. |5] 
Vvml(R) is somewhat different from the standard CML rotation 
curve. The difference is moderate for the "average" spiral as de- 
fined in this work, but in individual objects, it may not be so. 
From the wide range of colour profiles observed in individual 
galaxies (e.g. Figs. |2] and O, we expect cases in which the VML 
effect is more relevant than the average case considered here, 
while it will be negligible in other cases. 

Here, we take the function in Eq. |6] as the typical average 
contribution to the circular velocity of VML stellar discs. We 
combine it with a simple halo distribution, to mass-model the 
typical RCs of disc galaxies. Our main aim here is to compare 
the corresponding disc/halo decomposition to the results of the 
standard CML assumption; namely, we perform a differential 
analysis, rather than making absolute claims on disc fractions, 
halo profiles, etc. 

We consider the co-added rotation curves of Persic, 
Salucci & Stel (1996; hereinafter PSS96). These are ob- 
tained by binning a very large number of individual 
high surface-brightness, late type, bulge-free RCs (see 
also Salucci et al. 2007 for references on independent 
work on this issue. More information can be retrieved at 
http ://ww w.faceboo k.com/home.php#!/group.php?gid=3 102604] 
50630). One advantage of testing a template of co-added RCs 
with a template VML disc-velocity contribution is that this 
allows weak systematic VML kinematic features to be empha- 
sized also, which in individual RCs would be below the noise 
level created by modelling uncertainties and/or observational 
errors. On the other hand, relevant VML features occurring in 
a fraction of objects with strong colour gradients can only be 
studied in individual RCs. Ultimately, it is in the latter that the 
proof for the existence of VML discs must be found. 
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Fig. 7. Three representative examples of new VML disc rotation 
curves (solid and dashed lines) compared to a constant M*/L 
(CML) Freeman disc (dotted line). All velocities are normalized 
to the respective total disc mass and / band photometric scale- 
length. 



We assume the simplest possible dark-halo velocity contri- 
bution (see PSS96): 

V..«^ = Vi,(l-/^)(l.«^)^ ^ = ^ (7) 

with Rop, = 3.2 h[, Vopt the observed circular velocity at Ropt and 
(1 -yS)'^^ the fractional contribution to Vopi of the halo compo- 
nent. Correspondingly, the disc contribution (whether VML or 
CML) is normalized so that Vdis?(Rop,) = pvlp,. 

Although simple, the function in Eq. 7 has several mer- 
its. First, inside 4 h, it represents the NEW velocity profile for 
a < 0.2. Moreover, for appropriate choices of the parameter a, it 
also mimics the pseudo-isothermal and the Burkert (1995) halo 
velocity profiles. Then, it is worth recalling that, for any value of 
its structural parameters (a and /3), the slope of our halo velocity 
law (and of the profiles quoted above) is very different from that 
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Fig. 8. Disc/halo decomposition of the co-added rotation curves of PSS96 for different galaxy luminosities; the disc RC, the dark 
halo contributions, and the total RC are shown for each choice of VML (or CML) disc RC. Line symbols are as in Fig. |7j the dotted 
lines represents the PSS96 mass model with the CML disc curve. In panels (c,d) the total (disc+halo) RCs of the three VML disc 
decompositions overlap. 



of the disc contribution. At any radius we have 



d log Vl,alo(R; a,/3) d log Vdisc(R\P) 



dlogR 



dlogR 



> A 



dlogV(R) 



dlogR 



where the righthand side is the uncertainty with which we mea- 
sure the RC slope. This ensures that, in the case a) of high- 
quahty rotation curves (A\dlogV(R)/dlogR)\ < 0.04) and b) 
we know the actual distribution of the luminous matter, the 
mass model we obtain by x' fitting the RC is essentially unique. 
Values of a and disc mass significantly different from the best-fit 
ones, imply mass models that significantly fail to fit the observed 
rotation curve. In the mass modelling we neglect the gas contri- 
bution to the RC since it is always negligible out to ~ 3 h — the 
region where the differences between VML and CML RCs are 
maximal. 

We fit the co-added RCs Vcoadd(x) with our 2-free parameters 
velocity model: Vl^^^, = Vl,^(x; a,p) + Vl^^(x;/^), with Vvml 
given by Eq. |6] Let us recall that the co-added RCs we adopt 
here have been modelled in PSS96 with the same halo velocity 
profile of Eq. 7 plus a CML exponential thin disc. The mass 
models obtained there (shown in Figure 6 and in Eq. lla,b of 



PSS96) will serve as a gauge for the models we obtain with the 
present investigation. 

We perform analysis: the fitting parameters (a and /3) are 
optimized with a Leverberg-Marquardt algorithm (Press et al. 
1992). The increase in the r.m.s. scatter in the velocity data be- 
yond Ropi comes from the fact that, while out to R„p, the template 
velocities come from the co-addition of 616 individual curves, at 
outer radii they are obtained in a different, non-trivial way (see 
Fig. 3 and Section 4 of PSS96). Then, the related r.m.s. scatter 
is larger and consequently the velocity data for R > Rgpt do not 
contribute much toward setting the best-fit velocity model. 

The resulting best-fit solutions for 3 different VML cases, 
which progressively depart from the CML one, are given in 
Table 1 and in Fig.|8] compared to the CML solution. As a result 
of the true universality of the RCs of spiral galaxies, the r.m.s. 
errors of the co-added curves are very small: A VjV < 0.04 (for 
R < 3 h). Within the 2-component halo -i- disc model, this guar- 
antees that the Icr fitting uncertainties on the parameters of any 
velocity model that fits the data well will be also small. In fact, 
we find that the Icr contour fitting uncertainty on a and /? is be- 
tween 15% and 30%, not surprisingly similar to the correspond- 
ing uncertainties of the CML -i- halo mass model in PSS96. 
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Table 1. Parameters a and /3 and inner baryon dominance radius 
of the disc/halo decompositions in Fig. [8] for CML and VML 
disc RCs. 





< M, 


> = 


-19.4 


< M, 


> = - 


-20.9 




a 






a 






VML (s= 1.0) 


0.87 


0.17 


0.32 


0.92 


0.39 


0.76 


VML (s=0.7) 


0.80 


0.13 


0.25 


0.75 


0.30 


0.55 


VML (s=0.5) 


0.75 


0.10 


0.20 


0.61 


0.21 


0.36 


CML 


0.94 


0.24 


0.42 


1.22 


0.53 


1.08 




< M, 


> = 


-22.4 


< M, 


> = 


-23.1 




a 






a 






VML (s= 1.0) 


1.28 


0.64 


1.34 


1.29 


0.69 


1.53 


VML (s=0.7) 


0.94 


0.52 


1.05 


0.93 


0.57 


1.18 


VML (s=0.5) 


0.75 


0.42 


0.81 


0.75 


0.46 


0.91 


CML 


1.63 


0.80 


1.93 


1.71 


1.85 


>2 



First of all we notice that, when also considering the fitting 
uncertainties, the VML stellar disc contribution to the total RC is 
dilferent from the CML disc contribution. At its peak, the ampli- 
tude is lower by ~ 20%. Noticeably, the VML discs (in combina- 
tion with the PSS96 DM halo profiles) reproduce the kinematic 
data in an excellent way: the resulting values are not higher 
than those of the CML models in PSS96; that is to say, the VML 
hypothesis is by no means excluded by the data. With respect to 
the standard CML decomposition, the overall mass distribution 
of spiral galaxies inferred in the VML scenario is not drastically 
modified, but some major differences are present. These emerge 
beyond the uncertainties on the structural parameters, but rely of 
course of the fundamental assumption that Eq. 6 holds for most 
spiral galaxies. 

The VML disc mass distribution has a higher percentage of 
its total mass residing in the inner regions (i.e. 70% within R < 
2hf, Fig.|6l). The disc normalization, constrained by kinematic 
data, is lower and the inferred total disc mass and M*/L ratio are 
about 1.5 times lower than in the CML case. 

In low-luminosity spiral galaxies (Fig. [8^) dark matter al- 
ready dominates the mass distribution at small radii, so that fit- 
ting with VML discs does not influence the inferred halo proper- 
ties. At higher luminosity, however, the role of the disc becomes 
more prominent. Since the disc normalization is lower and the 
VML RC, after its peak atR~ l.l h, falls off more rapidly than 
the CML RC, in the region R > 0.5 Rop, the dark halo velocity 
contribution must rise more steeply, to match the (typically still 
rising) rotation curve. 

For galaxies of M/=-21 and brighter, the PSS96 mass mod- 
els showed a region of inner baryon dominance (IBD). Defining 
the corresponding radius Ribd as the region where the baryons 
contribute more than 50% of the observed circular velocity, the 
CML mass models yield Rjbd as wide as Rgp, or more (see 
Table 1). Mass models based on VML discs + PSS96 dark haloes 
also show an IBD region, but its size is smaller: Rjbd extends 
beyond Rop, only in objects of highest luminosity (Fig.[8j),c,d). 
In the luminosity bin of Milky Way-like galaxies (M/=-22.3, 
Flynn et al. 2006), the luminous and DM contribution to the 
rotation curve at Rop, are comparable, rather than still disc- 
dominated, as is the case for the CML mass models. 

We find, as in PSS96, that the parameters a, /3 and galaxy 
luminosity are closely related, so we expect that scaling laws re- 
lating disc and halo masses, luminosity, disc scalelengths, and 
other galaxy global quantities, still hold in the VML scenario. 
However, we do not attempt to give explicit fitting formulae 
(similar to those in PSS96) because the differences we find 



with the PSS96 results are moderate. Moreover, the VML sce- 
nario needs yet to be confirmed by combined kinematic/chemo- 
photometric analysis of individual objects — where scatter is 
large (Fig. |2] and [3]). 

For galaxies in the highest luminosity bin, there may be 
some weak feature in the co-added RC and in individual RCs 
that may be traced back to a VML disc. At 1-2 h;, the CML 
-Hhalo rotation curve slightly overpredicts V(R) at a 1 cr level in 
2-3 velocity data. The VML-i-halo model instead matches all the 
data (see Fig. 6 of PSS96 and Fig. [8};,d). Moreover, at R ~ 3h, 
VML-)-halo model RCs tend to remain flat, rather than decline 
like the CML-i-halo curves (Fig. [8};,d). This is not strongly ex- 
cluded by the co-added curves (due to the larger error bars be- 
yond Rop,) and could be in line with the RC profile of some fast 
rotators (Spekkens & Giovanelli 2006). 

The investigation of the core vs. cusp issue in spiral galax- 
ies requires the mass modelling of individual RCs; however it 
is worth discussing whether, with VML discs, the NFW haloes 
may become compatible with the observed kinematics. Mass 
modelling of co-added RCs with a CML disc + the "neutral" 
halo of Eq. 7 yields a core in the DM density profile; such core 
is similar to that found by the (much more decisive) analysis of 
the kinematics of selected individual objects. By replacing the 
CML disc with a VML one, we still find almost the same den- 
sity core; although its size is about 30% smaller, it is still large 
enough to exclude the NFW profile beyond the uncertainties in 
the fitting method. Thus, at face value, VML discs cannot recon- 
cile N-body predictions and observations. 

Furthermore, we can resort to the ingenious proof of incom- 
patibility with the observed kinematics given by Salucci (2001), 
to show that VML discs are more of a problem than a solution 
for NFW halos. The argument involves comparing the slope of 
the observed RCs with the values predicted by the NFW profile. 
We here apply this line of reasoning at the radius where the disc 
RC has its maximum, rather than at Rop, as in Salucci (2001). 
The following chain of facts and evidence leads to the inevitable 
conclusion that the DM distribution, at least in low luminosity 
spiral galaxies, is shallower than that predicted by a NFW halo: 

i) At the radius in which the disc contribution peaks (2.2 h for 

CML discs) and the disc RC has a flat slope, we have 

dlogVh , r, s-i^logV 

-7; —\2.2h - (1 - P2.2) — -\2.2h 

d log R d log R 

with 132.2 - iV^/V^)\2.2h the fractional disc contribution to 
at 2.2 h. 

ii) A generous estimate is ^62.2 ^ 0.15. 

iii) For NFW haloes it is dlogV^FW /dlogR < 0.2 for R > 
1.5/1. 

iv) For RCs of low luminosity objects (Mj < -21), we have 
dlog V/dlogR\2.2h > 0.3 (PSS96). 

In the VML mass modelling, the same chain of arguments also 
applies, with only the difference that now the peak radius is at 
inner physical radii, i.e. at ~ IJh. Noticeably, at this radius the 
observed RCs show an even steeper slope, dlog V/dlogR\i jh > 
dlog V/dlogR\2.2i„ and therefore the discrepancy is even larger 
than in the CML case. For instance, for objects of M/ ~ -21, the 
observed dlogV/dlogR\2.2 1, ~ 0.2: their RCs at 2.2 h can be, 
in the CML scenario, marginally reproduced by a NFW halo. 
But in the VML case where the disc contribution peaks, we have 
dlog V/dlogR\i j /, ^ 0.25, inconsistent with the NFW halo. 
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5. Summary and conclusions 

We have investigated the radial profile of the stellar mass-to-light 
ratio in spiral galaxies, and seen how their inclusion changes the 
picture of the mass models of spiral galaxies. We derived the ra- 
dial variations of M*/L in galactic discs, as predicted by chemo- 
photometric models in the "inside-out" formation scenario and 
as suggested by the observed colour gradients of spiral galaxies. 
The differential radial variation in M*./L we predict for a "typ- 
ical average" spiral turns out to be largely independent of the 
assumed SF history, stellar IMF, global galaxy colour, and other 
modelling details. It follows a profile of the type 



L, 



- exp < -a{s) 



with 0.5 < s < 1.1 and a(s) given in Eq. 3b. The predicted ra- 
dial M*/L variation can be significant, for instance in the / band, 
M*/L/ can decrease by up to a factor of two between 1 scale- 
length and the outskirts of the stellar disc (4/i). We provide the 
corresponding variable mass-to-light rotation curve Vvml{R), 
that describes our best guess of the stellar disc contribution to 
the circular velocity in a "typical average" spiral galaxy. Since 
small variations (~ 20%) in the latter quantity can also affect the 
mass modelling, VML effects should not be neglected in objects 
showing radial colour gradients. 

The resulting disc circular velocity Vvml{R) has some in- 
teresting features with respect to the classic constant mass-to- 
light (CML) Freeman curve. Since M*/L decreases outwards, 
Vvml{R) peaks at smaller radii: at 1.5-1.9 h; (rather than at 
2.2 h), and it approaches the Keplerian fall-off earlier, at 3-3.5 hj 
(rather than at > 4 h). For a given observed RC, the steeper de- 
cline of the VML stellar contribution at the edge of the optical 
disc, with respect to the CML case, implies there an even steeper 
radial rise in the halo contribution. 

With the average VML RC of stellar discs (combined with a 
cored dark matter halo), we have successfully fitted the co-added 
rotation curves of spiral galaxies of PSS96. The radial varia- 
tions in the stellar M*/L predicted by chemo-photometric mod- 
els seem compatible with presently available kinematics (Fig. [8]). 
Independent support may come from Salucci et al. (2008), who 
compared photometric and CML kinematic estimates of disc 
masses: despite the generally good agreement (the rm.s. and 
the offset being only of 0.2 dex and 0.1 dex respectively), there 
were a few cases of photometric masses higher than the kine- 
matic ones. These discrepancies could stem from VML discs. 

The VML disc scenario may be very relevant in two open 
issues of galaxy structure. Firstly, we have shown that the ra- 
dial variations of M*/L allowed by chemo-photometric argu- 
ments, cannot reconcile the DM halo density profiles derived by 
mass modelling with those emerging from ACDM simulations. 
In contrast, the "standard" assumption of constant M*/L helps 
NFW haloes to fit the observed kinematics. Because the colours 
of spiral galaxies typically redden inwards, the corresponding 
M*/L variations are an additional problem for NFW profiles. 

Secondly, we discuss whether, inside the optical disc, the 
stars alone, without any dark component, could account for the 
observed RCs via a radial variation of M*/L. This question was 
raised in the 80s, when evidence for DM on galactic scales 
started to gather. It is usually bypassed in view of the existence 
of non Keplerian, very extended HI RCs implying a huge mass 
discrepancy beyond the optical region. Nevertheless, the issue 
of the RC within the optical region is far from trivial, as no-DM 
scenarios like MOND and f(R) theories of gravity are becoming 
popular (e.g. Sanders & McGaugh 2002). 




Fig. 9. No DM, VML model of the co-added RC at M/ = - 19.4. 
The radial coordinate is in units of Ropt = 3.2 /?/. 



Figure |9] shows the co-added RC of galaxies with < M/ >= 
-19.4 out to 3.5 hj. It has a steep slope which implies a clear 
mass discrepancy from R = hj outwards. We recall that at these 
luminosities the HI disc is dynamically important only outside 
the region we are considering here (de Blok et al. 2008). In Fig.|9] 
we show that the co-added curve can be reproduced well without 
any DM, with just a VML stellar disc. To reach this, however, 
we must invoke an increase in M*/L with radius at the pace of 
1 + 0.35{R/hj)^. The corresponding colour variation (PST04) is 
a reddening of over 1 mag in {B - I) between R = 0.5 hi and 
3hj. This is totally at odds with the observed radial bluing of 
about 0.4 mag over the same radial range (Fig.[3]l. We conclude 
that a dark component must also be present inside the luminous 
regions of disc galaxies. 

The existence of VML discs may play an important role 
in the mass modelling within the MOND paradigm, where the 
observed distribution of the ordinary matter, stars, and HI gas 
is argued to completely account for the observed kinematics. 
In MOND, the circular velocity due to the distribution of the 
baryonic matter is "boosted" with respect to what we would 
have in Newtonian dynamics. From about 1 scalelength out- 
ward, a MOND self-gravitating Freeman disc implies a rotation 
curve constant with radius and with an amplitude that, at Ropr, is 
about 20% more than in the Newtonian case. We consider high- 
luminosity spiral galaxies with V„pi ~ 250 km/s, M/ ~ -23.5, 
/!/ ~ 5.5 kpc, belonging to the second highest velocity bin of the 
co-added curves in PSS96 and shown in Fig. [10] In these objects 
the dynamical contribution of the HI disc and the stellar bulge 
are negligible in the region between 2 and 3 disc scalelengths. 
For a VML disc the stellar surface density is given by Eq. 5. We 
model the RC of a VML disc in MOND dynamics with ^ = 1 
(shallow profile) or s - 0.5 (concave profile). Let us recall that 
in Newtonian gravity both VML discs, or a CML disc, can fit 
this same RC when combined with a suitable cored dark halo 
(Fig.[8}l). In MOND, instead, the outcome is different: while the 
i = 1 VML disc accounts for the RC, the s - 0.5 case does not 
(Fig. [TOb . This result is not a test of MOND, since we do not 
know the average s in spiral galaxies, but indicates that in alter- 
native theories to dark matter, the effects of M*/L variations are 
much greater and potentially testable. 

In the disc+dark halo mass modelling, the CML and VML 
discs yield equally good fits to the observed RCs (Fig.[8]l, so that 
one cannot use the RC to significantly reveal VML features. This 
is a non-trivial result: it is due not to an intrinsic degeneracy of 
the mass modelling, but to the fact that the allowed M*/L varia- 
tions are not very large. In the "average" late type spiral galaxy, 
chemo-photometric arguments indicate a decrease in M*/L/ by a 
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Fig. 10. MOND VML mass model of a co-added rotation curves 
of luminous objects, s = 1 (dashed) and s - 0.5 (solid) cases. 

factor of 2 over the optical region, which is hardly distinguished 
from the kinematics. In the simplest VML disc case (s = 1 in 
Eq. |5]), for instance, the surface density distribution remains ex- 
ponential, with just a scalelength about 20% shorter than that of 
the /-band luminosity profile. 

Noticeably, in alternative scenarios to DM where the kine- 
matics are completely determined by the distribution of the bary- 
onic matter, the mild radial M*/L gradients we predict are ex- 
pected to play a relevant role. 

Although the predicted M*/L gradients do not radically alter 
the main results of traditional mass modelling, our comparison 
of the VML and CML assumptions shows that VML discs have 
an impact on RC decomposition: the trend is to obtain smaller 
disc normalization and smaller, but non-negligible, dark halo 
cores. The NFW halo profile becomes even less compatible with 
the observed kinematics in the inner regions; therefore, M*/L 
gradients should in general not be neglected in mass modelling, 
and a combined chemo-photometric and dynamical analysis of 
individual objects (especially those showing strong colour gra- 
dients) is the best approach. 

More ambitious, but not beyond reach, is the idea of using 
detailed kinematic analysis as an independent test of the predic- 
tions of chemo-photometric models. By turning the above argu- 
ment around and by assuming a complete inner baryon domi- 
nance in the innermost regions, one can possibly aim at measur- 
ing the variations in the stellar M*/L directly from high-quality 
rotation curves and check the predictions of chemo-photometric 
models. Also, since the prediction is that the scalelength of the 
stellar density profile is shorter than the luminous scalelength 
(by at least 20%, 5=1 case), future investigations could aim at 
revealing VML discs by estimating the dynamical scalelengths 
from the observed RCs and then comparing them with the pho- 
tometric ones. 
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